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We study the one-loop effective action for gravity in a cosmological setup to determine possible 
cosmological effects of quantum corrections to Einstein theory. By considering the effect of the 
universal non-local terms in a toy model, we show that they can play an important role in the 
very early universe. We find that during inflation, the non-local terms are significant, leading to 
deviations from the standard inflationary expansion. 


I. INTRODUCTION 

Quantum corrections to Einstein action have received the attention of numerous authors for quite a long time. In 
the original works pj, the divergent structure of the theory was determined. It was seen that the theory required 
at the one-loop level OtJZ 2 ) counter terms (albeit it turned out to be finite on shell if matter was absent). These 
divergences will actually be important for our discussion since, by unitarity arguments, they determine the non-local 
part of the one-loop effective action. 

Although it was soon found out T] that higher order loops required more and more counter terms, making Einstein 
gravity non-renormalizable, not even on-shell, this does not mean that one-loop corrections are useless. In fact in 
a regime of small curvatures one-loop corrections will certainly dominate over two-loop corrections and so on. As 
a matter of fact this is not very different from the familiar expansion in chiral perturbation theory in powers of 
{p 2 /16nf 7T ) 2 (see e.g. 0 ). 

Within this philosophy, the one-loop effective action for Einstein gravity has received quite a lot of attention recently 
0. For instance, the quantum corrections to Newton law have been considered by a number of authors 0, and after 
some controversy the correct result has been found. Quantum mechanics gives a correction to the classical result of 
0(l/r 3 ) and positive; that is, at long distances gravity is more attractive that what Newton law predicts. Of course 
the correction vanishes anyhow at large values of r and it is accompanied by a very small coefficient, so for terrestrial 
or astronomical purposes this correction is most likely irrelevant and unlikely to be tested ever. 

However, the situation may be more promising from a cosmological view point. There is a cumulative effect of 
gravity and, given a fixed density of energy, the integration of this effect over large volumes could give an observable 
signal. This is, in short, the possibility that we would like to tentatively study in this work. 

We shall consider a flat Friedman-Robertson-Walker metric background and determine the effect of quantum cor¬ 
rections on the cosmological evolution of the scale factor. We shall concentrate here on the consequences of the 
non-local terms that necessarily appear in the effective action due to unitarity considerations. It is a well known fact 
from pion physics that these terms dominate al large distances, but somewhat surprisingly they do not appear to have 
been considered before in the present context (with one exception to which we shall turn below). Somewhat to our 
surprise, the effect of these quantum corrections turns out to be relevant. 

In the present, exploratory, paper we shall not consider the full one-loop effective action, but shall limit ourselves to 
a toy model where only the scalar curvature is included in the effective action and neglect the full Ricci or Riemann 
tensors. We shall also make for simplicity a number of additional approximations that we shall discuss in more detail 
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below. The action of our toy model up to TZ 2 terms is 


S = j dxy/^g (k 2 K + alZ ln(V 2 //j, 2 )7?. + [31Z 2 ) + S M 


where k 2 = 1/(167tG) = Mj> ; /(167r) and Sm includes the matter fields (and in particular the inflaton sector of the 
theory). Note that the constants a and (3 are dimensionless. The expansion parameter is curvature , TZ, which in the 
approximation H = 0 is in a simple relation to the Hubble parameter, TZ ~ H 2 . Hence, when H/Mp\ <C 1 e.g. at 
present times, local higher order terms can be neglected. The values of a and f3 are on a different theoretical footing. 
a is entirely determined from the analytical structure of quantum corrections induced by the lowest dimensional term; 
in other words, its value is uniquely determined once one insists in the long distance description of gravity being 
provided by Einstein theory. /?, on the other hand, is model dependent. Its value is fixed as a boundary condition 
upon integration of other degrees of freedom that do not have been included in ©. Furthermore, (3 is scale dependent 
so as to cancel the (local) In /i 2 dependence appearing normalizing the log in (JTJ . 

We neglect higher orders in the expansion in derivatives of the metric. A power counting can be established here 
5), in parallel to what is done in chiral perturbation theory. We shall therefore work with a precision where only up 
to four derivatives of the metric need to be included. 


II. DETERMINATION OF THE EQUATIONS OF MOTION 

We split the action (ignoring the matter part for now) into three parts and redefine the constants for convenience 

S = „ 2 (/ dx^/^gTZ + a J dxyf^g TZln(V 2 //j?)TZ + {3 J dx\f^g TZ 2 ^ 

= k 2 (s 3 + &S 2 + PS^, (2) 

where a = n~ 2 a, (3 = n~ 2 f3. The dimensionful constant ^ is actually a subtraction scale that is required for 
dimensional consistency. The coupling /3 is p. dependent in such a way that the total action S is ^-independent. In 
conformal time, dt = a dr , we have 


9pv = a 2 {T)ri^, 1Z = -6^, V 3 9 = aV)- 

a J (r) 

Variation of the local action is straightforward and gives the following results 
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The variation of £>2 is more delicate and requires some discussion. The d’Alembertian in conformal space is related 
to the Minkowski space operator by @ 


1 


V 2 = a + -K. 
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To the precision we are working in the curvature expansion we can neglect the TAterm in the expansion of the 
d’Alembertian and commute the scale factor a with the flat d’Alembertian; therefore we set 


V 2 = (-) □, 

Va 0 / 


( 7 ) 


Where ao = a(0). The rescaling (absorbable in (3) ensures that at r = 0 the d’Alembertian matches with the 
Minkowskian one. In fact, from now on we will set ao = 1 for simplicity. 

We can now separate S 2 into a local and a non-local piece 


S 2 = j dx\f^g (^~2TZ\n{a)TZ + TZhi{U\/gL 2 )v}j 


= S^+S. 
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Variating the first term is again straightforward 
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We shall discuss the variation of S^ 1 next. 


A. Evaluation of the non-local contribution 


The genuinely non-local piece in S 2 is 


S. 
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J dx\/~^g lZln(0/y 2 )lZ 

J dx\f-g{x)ll(x) J dy\/—g(y) (z| \n{\3/y 2 )\y)K(y). 


(9) 


To evaluate this logarithmic term, we make use of the identity ln(ir) ss — l/e + a; c /e, valid when e is small. Using this, 
we see that we need to compute e~ 1 (x\{D/y 2 ) e \y}, that has the integral representation 


' 1 (a, , |(D/^ 2 ) e |y) = -2n 2 y 2e / dk k 2+2e 
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where J\ is the Bessel function. As we are only interested in the time evolution of the scale factor we can integrate 
out the spatial dependence in Eq. absorbing the constants into d leaving 


s¥ = J dTy/-g(T)K(T) J dr' y/—g(T') y 2e _ (11) 

The above integral is of course divergent; the divergence is local and is exactly cancelled by the — 1/e in ln(a:) ~ 
— 1/e + x e /e. By taking this term into account and taking e —> 0. © can be easily numerically computed for a given 
background. 

A final word of caution has to do with the causal conditions to be imposed on our Green functions. We have been 
so far rather cavalier about that. When deriving evolution equations (as opposed to S'-matrix elements) one has to be 
careful with the causal conditions. The usual Feynman rules lead to the so-called in-out effective action, relevant for 
S'-matrix elements. When interested in causal evolution one has to consider the in-in effective action This affects 
the causal definition of the Green function A(s, y ) = (x\(D/y 2 y\y}. If Ap, Ad, A + and A_ are the Feynman, Dyson, 
advanced and retarded Green functions, respectively, the combination that actually appears in the in-in effective 
action takes the matrix form (see e.g. 0 for a clear discussion on this formalism) 


(n + n.) 



( 12 ) 


where the subscripts denote positive and negative frequency, respectively. This ensures that when taking a functional 
derivative w.r.t. a(t), only earlier times contribute to the evolution. This can be enforced by setting the limits in the 
time integrals appropriately. 


B. Variation of the non-local contribution 

Now we can proceed with the variation of the non-local part (for simplicity, we set e = 0 here; although the 
expressions are all ill-defined, we can replace e back easily) 




















The derivatives can be easily performed by integration by parts 



+ J dr' —— -{a"(r ') 2 + 2a , (r , )a"'(r / ) + a(r')a (4) (r')) . (14) 


The poles at r = 0 terms are an artifact arising from the fact that at r = 0 we patch together Minkowski space and 
de-Sitter space by starting inflation at that point. If this is done smoothly enough, the derivatives of the scale factor 
vanish at that point (this is evidenced by the fact that all these terms contain derivatives of a at r = 0). Hence, we 
disregard these terms in the following. We have checked that modifying the matching has unobservable consequences 
in our results. 


III. EFFECTS ON INFLATION 


With the variated action, we can now look for solutions of the resulting equation. In general, the equation is a 
complicated integro-differential equation that is most suited for numerical studies. However, we expect that the new 
terms may have an effect during inflation when TZ is large and H is small. 

In conformal time, the scale factor grows during inflation as 


a/(r) 


1 

1-Ht 


(15) 


where H is the inflationary Hubble rate as in a = ao exp (Ht). It is evident that in conformal time we must restrict the 
time interval considered to r < l/H. Note that in order for this to be the solution of the varied gravitational action 
(without the higher curvature terms) one has to put in the appropriate matter terms into the action. Assuming a 
simple generic inflationary model, one can relate the the Hubble rate to the inflaton potential by H 2 = 8nGV((f>)/3 , 
where <j> is the value of the inflaton during inflation (approximated with a constant). Here we consider H to be a 
constant and study how the new terms affect the inflationary expansion. In a more realistic treatment left for future 
work, one has to let the inflaton roll and solve the coupled system of equations of the inflaton equation along with 
the modified Einstein’s equation. 

We proceed by solving the varied gravitational action by a perturbative approximation, i.e. we consider the non¬ 
standard terms as a correction to the standard inflationary solution. This procedure is only valid as long as the 
correction is small compared to the unperturbed solutions, which sets the limits for the validity of our approach. 

For calculations it is useful to note the relation 


a'j(r) = Haj(r), 


which implies that 

4 n) =n\H n a n I + 1 (r). 

The 0-th order equation, corresponding to including the inflaton sector in the action, is 

- 12a" + 24 H 2 a 3 = 0, 


which Eq. is a solution of. 
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A. Perturbing the 0-th order solution 


By using the 0-th order solution we see that the variation of S 3 vanishes (this is not a general result, i.e. this 
only happens with the 1Z 2 term, other terms of the form 7 Z n , n > 2 are non-zero). The variation of S 2 does not vanish 
and after a straightforward substitution we get for each part 
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Hence, the 1st order equation of motion for a(r) is 


( 20 ) 


a" — 2H 2 a 3 = 125d -8 H^a 


{ 

+dH 4 a 3 j [ dT'^^ + 8H 4 aj [ dr'+ 20H 4 ai [ j. (21) 

J o T ^ r 7o t - t' J o t - r' J 


The r.h.s. of this equation (that depends on the unperturbed 0-th order solution only) can be regarded as a driving 
external force for the evolution equation of the conformal factor. 

The integral terms can be explicitly computed (we now restore fi and e) 


dr 1 
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where F is the hypergeometric function. In the limit e —> 0, Eq. iPl) tends to 


dr' 


(jT^ (^ | 

( t _ 7 r ^i +2e = —-E(n, 1,1, Ht) + ln(Tfi)F(n, 1,1, Hr) + d e F(n, 1,1 — e, Hr) 


e—>0 


=-(1 - Hr)~ n + (1 - Hr)~ n In(r/x) + <9 e F(n, 1,1 - e, Hr) 
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e—>0 


(23) 


indicating the need to regulate the integrals by the addition of an appropriate counter-term (this is the role of the 
— 1/e term appearing in the logarithm representation). We are now ready to solve Eq. 11211 1 numerically, using the 
0-th order solution m as an initial condition. 


B. Numerical analysis 


As we can see from the form of the inflationary expansion in conformal time, the scale factor is singular at r = 1 /H. 
Furthermore, if we wish that e.g. ai/ao ~ e 60 ~ 10 26 , we must require that 1 — Ht ~ 10~ 26 , i.e. that r is very close 
to the singular point. It is hence useful for numerical work to do a change of variables, e s = 1/(1 — Ht). In these 
coordinates the inflationary expansion is simply a/(s) = e s and the coordinate s £ [0, oo]. The equation of motion 
with the new dimensionless time coordinate is (after regularization) 

e 2s a" + e 2 V - 2a 3 = 12 aH 2 ( - 8e 3s + 4e 3s G(4, s) + 8e 2s G(5, s) + 20e s G(6, s)) (24) 


where we have defined G(n, s) = e ns ln(-^(l — e s )) + d e F(n, 1,1 — e, 1 — e s ) 


and divided both sides by H 2 . In 

e —>0 


this form of the equation of motion, we note a number of interesting properties. First of all we see that the coupling 
constant appears in a combination aH 2 , indicating that the corrections to the standard evolution become extremely 
small at late times (recall that presently H ~ 10~ 42 GeV and that a is dimensionful and proportional to the inverse 
Planck mass). Secondly, the arbitrariness in the choice of the scale is exhibited by the presence of the n/H term inside 
the logarithm. Finally, we note that the right hand side is singular at s = 0 (r = 0) but this is only a logarithmic 
divergence and can be avoided in numerical work by starting the calculation at some small non-zero value of s. 

For numerical work, we need to estimate the scale of inflation, V(<f>i n f) = Vo, and hence the Hubble parameter 
during inflation, H 2 = 8nGVo/3. As an absolute lower limit, the energy density in the inflaton at the end of inflation 
must be enough to reheat the universe to a high enough temperature for nucleosynthesis to occur (T ~ 1 MeV), but 
typical values can be much bigger than this, up to the CMB normalization limit < 10 16 GeV, corresponding to 


H ~ 10 13 GeV. 


It should be stressed here that the model proposed in 0 is not a realistic one inasmuch as we only include IF 2 
terms and neglect other possible curvature contributions. Our purpose here is to test in a simple setting whether 
these quantum contributions could be important at all. Therefore, albeit the actual values of a are known (see e.g. 
0 ), after integration of gravitons and the rest of massless particles in the Standard Model (massive particles are 
irrelevant for this discussion as they do not provide logs), we provide numerical results for several values of a around 















(a) /x = Mpi. 


(b) /x = 1 meV. 


FIG. 1: The scale factor relative to the inflationary expansion for different values of a, /x and Vo. 


the value |ct| ~ 10~ 4 , which is the natural order of magnitude expected. A more accurate treatment would require 
the introduction of all the Riemann curvature tensor components. It would also require to establish which particles 
are exactly massless (or have a mass much smaller than the inverse horizon radius for that matter). 

There is also a built-in /x dependence in our results. As we have stressed, the results should be, in principle, 
/x-independent (this is, incidentally, a point that is often overlooked in this type of analysis). The dependence on // 
in the non-local piece is exactly compensated by the (logarithmic) /x dependence of (3. Of course we do not know 
the value of f3 as it contains contributions from all modes that have been integrated out and that are not explicitly 
included in the lagrangian. We can, however, estimate the ‘natural’ scale for /x; i.e. the one that minimizes higher 
order corrections. 

Taking into account that massive particles are integrated out and that they do not generate logs, within a renor¬ 
malization group approach it is natural to consider that the natural scale is that of the lightest particle that has 
been integrated out. This is similar to what is done in effective lagrangians for the strong interactions such as chiral 
lagrangians, where the optimal scale is somehow related to the scale of chiral symmetry breaking that separates ‘light’ 
degrees of freedom from ‘heavy’ ones. In the Standard model we assume this lightest mass to be 1 meV, corresponding 
to the neutrino. 

To understand the effect of all of the previous choices, we have solved Eq. (EH) numerically for different values of 
a , Vo and /i. They are shown in Figs^a) andQJb) as a ratio of the scale factor to the inflationary expansion, a/a/. 

The normalization scales are chosen to represent the minimum and maximum values. Note that here we only show 
results for positive a. For negative alpha, the curves simply turn in the opposite direction at the same value of s and 
therefore we choose not to shown them here. 

From the figures one can see how a higher inflationary scale leads to deviations from the inflationary expansion 
earlier than a lower scale. This is as expected since the source term in Eq. E3 is proportional to H and hence to 
Vo. Similarly, a larger a has the same effect. The effect of changing the normalization scale /i has a mixed effect. 
If // = Mpi , the ln(/x/i?)-term in G(s,n) is positive for all considered values of Vo- However, if /x = ImeV, the 
logarithmic term changes sign from negative to positive at Vo ~ 2 x 10 3 GeV so that at large Vo, the source term 
is negative. In the source term, the —8e 3s term is subdominant compared to the other terms. The hypergeometric 
function gives a negative contribution for all n values considered here. Hence, the the sign of the logarithmic term is 
crucial, as is clear from the figures. 

The relative contributions to the action arising from the logarithmic term is shown for a particular case in Fig. EJ 
along with the corresponding numerical solution. From the figure one can see that the perturbative approximation 
we are making is appropriate as the action is still dominated by the Einstein term, 1Z. Note that at then end of the 
corresponding calculation the relative scale factor has decreased to about 0.96. The new terms are also subdominant 
for the other choice of parameter values shown in Fig. ^ 
































FIG. 2: The scale factor relative to a/ (solid line), log 10 (|27?.z h\(a)lZx/lZ\) (dotted line), log 10 (|7?.i ln(D//r 2 )7?.i/7?,|) (dashed 
line) for a = 10 —3 , fi = ImeV and Vg= 10 16 GeV. 


IV. CONCLUSIONS 

In this paper we have considered the effects of quantum corrections to gravity in a somewhat simplified setting. 
By considering the non-local terms that necessarily appear due to unitarity considerations in any effective lagrangian 
involving massless particles, that typically dominate long-distance physics, we have found that during inflation the 
non-local effects are important and lead to deviations from the standard inflationary expansion. The effect is sizeable, 
as for typical inflationary parameter values the expansion rate is changed after only 10 — 20 e-foldings. The sign of 
the effect depends on the parameter values, and in particular on the sign of a and size of the normalization scale 
fi. Recalling that a change in /i is tantamount to a change in the coefficient of the local 1Z 1 2 3 4 term, pinning down the 
physical value for f3 (and its proper /.t dependence) is very important. 

Taking into account that quantum corrections actually strengthen gravity at long distances, we believe that in the 
physically relevant situation, inflation would be slowed down or halted by the quantum corrections. 

This type of effects have not, to our knowledge been considered before, except for the studies presented in |T7il j. 
Although no detailed numerics are presented in this reference, the authors conclude that quantum effects slow inflation. 
Unfortunately the two approaches could hardly be more different and hence comparison is hopelessly difficult. It would 
be nice to make a clear contact between the two approaches. 

In doing the calculations, we have done a number approximations and simplifications in order to see whether the 
quantum effects can be important. As this has proven to be so in the toy model considered here, the effects of the 
non-local terms need to be studied more carefully in a more realistic model. We certainly believe that this issue 
deserves further studies. 
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